Effects of Interlayer Interaction on the Superconducting State in YBCO 
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For a two layer system in a weak coupling BCS formalism any interlayer interaction, regardless 
of its sign, enhances the critical temperature. The sign has an effect upon the relative phase of 
the order parameter in each of the two planes but not upon its magnitude. When one of the 
planes has a dispersion consistent with CuO chains and no intrinsic pairing interaction there is both 
an enhancement of the critical temperature and an s + d mixing in both layers as the interlayer 
interaction is increased. The magnetic penetration depth, c-axis Josephson tunneling, density of 
states and Knight shift are calculated for several sets of model parameters. 
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I. INTRODUCTION 

The search for the mechanism which causes supercon- 
ductivity in the copper oxide materials is an ongoing ef- 
fort which has yet to reach a consensus. One factor which 
any model should account for is that the critical temper- 
ature tends to be higher in systems with multiple adja- 
cent Cu02 layers; and even in systems, such as YBCO, 
in which a CuO layer is adjacent to the Cu02 layers Tc 
seems to be enhanced. 

It is generally believed that the superconducting con- 
densate resides in the Cu02 planes, although one in- 
terpretation of the observed large x-y anisotropy of the 
zero temperature magnetic penetration depth (a factor 
of ~ 1.6) in YBCO indicates that there is as much con- 
densate in the CuO chains as in the Cu02 planes (ie, the 
condensate in the chains only contributes to the pene- 
tration depth for current in the direction parallel to the 
chains, ie the 6-direction) . [|l| Since it is believed that 
whatever mechanism is responsible for superconductiv- 
ity in the copper oxides is intrinsic to the Cu02 planes 
some other mechanism for creating superconducting con- 
densate on the CuO chains is required. 

In this paper we derive a Hamiltonian for a layered 
system and, making a simplifying assumption that there 
is no pairing between electrons which reside in different 
layers, derive a pair of coupled BCS equations for a system 
of two layers, each with possibly different dispersion and 
pairing interactions. In this model Cooper pairs can scat- 
ter between the layers so that, as in a two band model, 
even if there is no pairing interaction in one of the 
layers there will still be a condensate in that layer due to 
the interlayer interaction. 

Although we make a particular choice for the pairing 
interaction (which is motivated by the nearly antiferro- 
magnetic Fermi liquid model which leads naturally to a 
d-wave gap for single Cu02 planes. ^) a d-wave solu- 
tion is also found for other types of non-isotropic pair- 
ing interactions. One result of having gap nodes cross 
the Fermi surface is a low temperature linear behaviour 
of the magnetic penetration depth, A^^, as is observed 
in YBCO. jl|,^-0 Here we only try to model the low 
temperature behaviour and relative magnitude ||l| in the 



X and y directions of the magnetic penetration depth. 
The T ~ Tc behaviour is only reproduced for values of 
2Z\max/Tc about 1.5 times higher |^ than the value of 
4.4 found in the BCS weak coupling approximation. In 
the two layer model that we study here a higher value 
of 2Araa,x/Tc is obtained which is closer to that observed 
in experiments such as ARPES |^ or current-imaging tun- 
neling spectroscopy (ciTS) that measure the absolute 
magnitude of the order parameter. 

We find that the presence of the chains destroys the 
tetragonal symmetry of the Cu02 planes and shifts the 
d-wave gap nodes in the Cu02 plane off the diagonals 
in agreement with an earlier model, pd] In this case 
the gap contains an admixture of s- and d-wave sym- 
metry. Calculations of the c-axis Josephson tunneling 
current show that the positive and negative parts of the 
order parameter do not cancel, as for d-wave pairing in 
a tetragonal system, and that the Josephson junction 
resistance-tunneling current products, RJ(T = 0), are 
in the range of 0.1-3.0meV, in agreement with the ex- 
periments of Sun et al, [Q. We also calculate both the 
normal and superconducting density of states (dos) for 
the CuO 2 planes and CuO chains separately since some 
surface probes, such as CITS, |l^ can measure them sep- 
arately. Finally the Knight shift is calculated separately 
for both the planes and the chains. 

In section II we introduce our Hamiltonian and derive 
a set of coupled BCS equations for planes and chains and 
other necessary formulas, particularly the expression for 
the magnetic penetration depth in this model. In section 
III we present the the solutions of these BCS equations as 
well as the results of calculations of the magnetic penetra- 
tion depth, electronic density of states, c-axis Josephson 
tunneling, and Knight shift. Section IV contains a short 
discussion and conclusion. 



II. FORMALISM 

In this section we will present a Hamiltonian in which 
multiple layers are coupled through the pairing interac- 
tions between adjacent layers. We will then make the 
assumption that there is no interlayer pairs (ie, that 



1 



each Cooper pair resides in only one of the layers) and 
that there is no single particle interlayer hopping. The 
Hamiltonian of this special case is then diagonalized and 
coupled BCS equations derived. We then give expres- 
sions for the magnetic penetration depth in this model, 
the Knight shift, the quasi-particle density of states and 
the DC Josephson junction resistance-tunneling current 
product for a c-axis tunnel junction. 
The general Hamiltonian is: 



H 



,"/3 (akT,Q^kT,/3 + al^ f^aki^a^ 
^ X! ^k,q,a/375ak|,QaLki,/3«qT,7"-qi,5' (1) 
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where the greek indices enumerate the layers, the aj^ ^ 
(flk.a) create (destroy) electrons of momentum k in layer 
a (k is in units of where a is the lattice parame- 
ter), ek,a/3 is the electron dispersion, and Vk.q. 0/3-^5 is the 
pairing interaction. 

Performing a mean field analysis, we get: 

H = ^ ek,Q/3 (akt,Q"kT,/3 + akx,/3"ki,Q 

k,a/3 

~ (^k,a/3a|,|,„aLki,,3 + H.C, 

k,Q/3 



where C is a constant, H.c. indicates the Hermitian con- 
jugate, Z\k,a/3 = H;,q,a/375Xq,75 are the order 

parameters and Xq,7<5 = ('^^qT.i'^-qi,'?) ^'^^ "P^^^ 
ceptibilities. 

Writing this in Nambu spinor notation, we get: 



k,Q/3 



where At 

K.I 



'^kT.Q '''-ki^Q 



and: 



£W,a0 ^k,a/3 



(2) 



In this model the magnetic penetration depth is given 
by the expression: |13 



±i5k,a are the eigenvalues of Eq. (g), f(x) is the Fermi 
function, 5ap is a Kronecker delta and Uk^af3 is the uni- 
tary matrix which diagonalizes Eq. (^. The quantities e, 
h, and c are the electron charge, Planck's constant and 
the speed of light respectively. 

We now make the following simplification: T4;.q.Q7 = 

^^,q,Qa77 — ^^,q,Q/37(5f^a/3^7(5 and £^k,a = ^k,aQ — 

Sk,af3Sa[3- This meaus that there is only intralayer pair- 
ing and no interlayer pairing (ie, Z\k,Q = '^k,a/3Saf3 and 
Xq^^ = Xci,'ysSaf3 are both diagonal in the greek indices) 
and there is no single particle interlayer hopping. This 
Hamiltonian has the same form as that for a two band 
model studied by one of us |^ in an earlier publication 
and is similar to that studied by others. ]l^-p^ Interlayer 
pairing [|6|-0 has also been studied. 

The Hamiltonian has eigenvalues given by i?k,Q = 

±^e^ a + ^k a diagonalized by the unitary ma- 

trix: 



C^k,a 



"k,a Wk,, 



(4) 



where: 





Ek.i 



-Bk,a 



are the usual BCS coherence factors. Using this unitary 
transformation we can evaluate the pair susceptibili- 
ties to get: 



Xq,ct — \1qt,a^— qi,Q/ 

= ^fL^tanhf^^^ 



(5) 



where T is the temperature and is Boltzmann's con- 
stant. Note that if we had included the interlayer pair- 
ing from Eq. ^ we would have susceptibilities of the form 
(aqT,Qa-qi,/3) with a ^ (3 and both the eigenvalues and 
the unitary matrix would be much more complicated. 

For a bilayer system (ie, a = 1,2) the BCS equations 
are: 



E T-k.aA'^a [G^Ma^o - ^k,a/3j (3) = ^Y] (M.,q,llXq,l + M.,q,12Xq,2) 



where: 



G 



k,Q/3 



7k,Q/3 



9£^k.Q 

/(£'k,a) - f{E\^,p) ^-y _ ^ 

7(5 



^k,2 — ^ X] (^k,q,12Xq,l + ^k,q,22Xq,2 ) , 



(6) 



where we have taken Vk,q,i2 = Vk,q,2i, although in gen- 
eral only Vk_q_i2 = 21 is required. 

Noting that Xq^2 changes sign (see Eq. ||) with Z\k,2 
we see that this set of equations (g) is unchanged by the 
substitution {/\k,2, Vk,q,i2} {-^k,2, -Vk^q^ia} which 
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means that the overall sign of Vk_q42 only affects the rela- 
tive sign of the order parameters in the two layers and not 
their magnitudes. This is interesting because it means 
that the effect on Tc of having an interlayer interaction 
is independent of whether this interaction is attractive 
or repulsive, although some calculated properties (eg, c- 
axis Josephson tunneling current) still depend upon the 
relative sign of the interlayer interaction. It is important 
to emphasize that any interlayer interaction, either at- 
tractive or repulsive, tends to enhance Tc and that this is 
consistent with the observation that is higher in mate- 
rials with multiple adjacent CuO layers. This well known 
result can be easily shown by examining the coupled BCS 
equations ^ near T ^ Tc- In this limit we can write: 

where and Vq,^ are numbers and 77^ is a normalized 
function which could be taken to be d-wave and corre- 
sponds to the highest T^. The coupled BCS equations (6) 
can then be written as: 



iiVn^V^tanh 



£q,l 

2M 



n ^ En. 2 V 2fcR 



A. = 



gq,2 

2fcBT, 

£q,2 



(7a) 



l^i^tanh 

U ^ £q^2 \2kBlc 



(7b) 



We now assume that A^ is the dominant superconduct- 
ing channel when V12 — V21 — and obtain assuming an 
infinite band with cutoff ujc'- 

A = (AA,,+Z\.A,,)ln(^i:^^ (8a) 

A, = {A,X,, + A^K,) In (^^^^ , (8b) 

where Ay = Vij x the density of electronic states at the 
Fermi surface. Substitution of equation (8b) into ( |8a| ) 
leads to a quadratic in In {I.13ujc/Tc) with solution: 



Tr = 1.13e^/^ 



with: 



A = 



All + A22 + V (All — A22)^ + 4A12A21 



(9) 



(10) 



This result is well known and is given by equation (6.3) 
on page 105 of Mechanisms of Conventional and High Tc 
Superconductivity. |23| It is also found as equation (40) 



of H. Chi and Carbotte. g] We note that A12, whatever 
its sign, increases A and so increases Tc- 

If we had taken T4;,q,i2 as complex the symmetry 
would be {Zik,2, Vi,q,i.} ^ {/ik,2e-"^,yk,q,i2e"^} where 
^k,q,i2 = |Vk,q,i2|e*''^, and the relative phase between the 
layers would no longer be ±1. 

We note that by performing the unitary transformation 
S^HS where: 
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and making the substitutions Vu — V22 = + 
V12 — — V±, El = e + t and £2 = e — t we obtain 
both the Hamiltonian and BCS equations used by Liu et 
al. p9[ | Our work differs from theirs in that we allow both 
the dispersion and the interaction to be different in the 
two layers. This is important not only because we are 
able to model systems such as YBCO in which there are 
Cu02 planes and CuO chains, but also because the order 
parameter in each of the layers may differ in sign even 
in two identical layers. Other workers have studied 
models in which the electrons in the pairs reside in dif- 
ferent layers (ie, in which only Xk,i2 is non-zero) as 
well as models in which there is no intralayer interaction 
Jl6| , p^ (ic, in which only Vk,q,i2 is non-zero). 

After solving the set of coupled BCS equations (|^) at 
T = using a fft technique |^,|ll],|oj we approximate 
the order parameters, Z\k,i and Z\k,2, with: 

A^,c.^[A'^:^ht^ + At\t^ + A^\\^^) 

X tanh [l.lA^jTc/T- l) , (11) 

where the 77^^ are the three lowest harmonics given by: 



(So 

'7k 



= 1 



'7k = cos(A;:r) + cos{ky) 
Vu^ = cos(A:j,) - cos(fcy). 



and the Aa^ are their amplitudes. The amplitudes of the 
higher harmonics are all very much smaller in magnitude 
and the gap nodes and the maximum magnitude of the 
gap, which are the most important features of the order 
parameter, are essentially unchanged by this approxima- 
tion. We also calculate the magnetic penetration depth 
which in this system, since the 7^'|j^ are diagonal in the 
greek indices, is given by the simplified expression: 

2 _ 47re2 1 ^ a^k.Q de-k,a 



dki dkn 



k,Q 

/a/(£k,a) ^/(-Bk.a) 



(12) 
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which is the usual expression pO[ summed over the layers. 

The curvature of the penetration depth curve, A^^(T), 
(and also its low temperature slope) is governed by the 
ratio 2Z\max/7c, where Z\max is the maximum value of 
the order parameter in the first Brillouin zone, and 
is close to a straight line for the d-wave BCS value of 
2Z\max/7c =4.4. The presence of the interlayer interac- 
tion increases this ratio and makes the A^^(T) curve have 
a downward curvature. Experimental measurements of 
both the ratio 2Z\,„ax/7c |lO| , pT| as well as the penetra- 
tion depth in high quality crystals of both YBCO j^] and 
BSCO indicate that this ratio is quite high in the 
HTC materials - on the order of 6 or 7. 

Other quantities calculated are the Knight shift which 
is given by: 



i^(T)oclE 



dEu ' 



(13) 



the normal state electronic DOS which is given by: 

k 



lim 



1 



r 



r^o nil ^ (£k - + ' 



(14) 



the superconducting electronic DOS which is given by: 
1 



k 



lim 



1 



r 



r^o nn ^ (Ek - LuY + 



(15) 



and the c-axis Josephson junction resistance-tunneling 
current product, RJ{T), through a superconductor- 
insulator-superconductor junction for incoherent c-axis 
tunneling is given by the relation [p4|: 



RJ{T) = 
where: 



2t:T 



,L(R) 



in which the superscript L(R) indicates on which side 
of the junction the dispersion and order parameter are 
on, the sum over Un = 7rT(2ri — 1) is for all Matsub- 
ara frequencies, R is the resistance of the junction and 
7V^(^)(0) is the normal state electronic DOS given by 
equation (^^. If the tunneling were coherent the ma- 
trix element (which is incorporated into R) would have a 
(k — k') dependence, and the sums over k-space wouldn't 
be separable. 



III. RESULTS 

In this section we make an explicit choice for the dis- 
persions and interactions and then present the results of 
our numerical solutions to the coupled BCS equations (||) 
as well as the results of our calculations of the magnetic 
penetration depth, densities of states. Knight shift and 
Josephson current. As we wish to model YBCO we will 
want to account for both the Cu02 planes as well as the 
CuO chains. Further, we will assume that we do not have 
a pairing interaction in the chains, but only in the planes 
as well as an interlayer interaction. This means that all 
of the order parameter in the chains is due to the in- 
terlayer interaction. We note that although our solution 
technique [Q allows the order parameters to be complex 
and to have a relative phase between layers we find that 
in the models studied here, to within an overall phase, 
the order parameters are all real with a relative phase of 
±1. 

For the dispersions, Ek.a, we use: 

£k,a = - 2ta [(1 + Eq) C0s(fc2;) + COs{ky) 

- 2Ba coa{k^) cos{ky) - (2 - 2B„ - A*a)] , (17) 

where the parameters {^0,60, -Bq, /ia} are chosen so that 
the Fermi surface and bandwidth are close approxima- 
tions to those observed experimentally |^ (see Fig. |^). 
In order to model YBCO we chose {100, 0, 0.45, 0.51} for 
the planes and {—50, —0.9, 0, 1.2} for the chains. In both 
dispersions ta, which sets the overall energy scale, is in 
units of meV. For the interactions, Vk_q_ct^, we chose an 
MMP form: [01 



-Xo 



l + ^2|k_q_Q|: 



where Xo is a constant that sets the scale of the suscep- 
tibility, is the magnetic coherence length, Q = (tt, tt) 
is the commensurate nesting vector, and gafs is the cou- 
pling to the conduction electrons, the size of which can 
be fixed to get a desired value of the critical temperature 
and can be considered to contain Xo ■ The remaining pa- 
rameter, ^o, is given in reference j2| and will not be varied 
in this work. In this paper we set §22 = ie, there is no 
intrinsic pairing in the chains. This means that any su- 
perconductivity in the chains is induced by the interlayer 
interaction, 1712, since we have set the hopping between 
layers to zero. The effect of an interlayer hopping has 
been extensively studied in works by Atkinson and Car- 
botte |l|l as well as others. |l|Jl|,|l3-|9| 

We solve the coupled bcs equations (^ using a fft 
technique. In Fig. || we plot the lowest three Fourier 
components ( |ll[ ) of the zero temperature order parameter 
as a function of the interlayer interaction (higher Fourier 
components are all much smaller in magnitude) for two 
different intralayer interactions (upper and lower frames) 
in the planes (left frames) and chains (right frames). The 
values plotted in Fig. || are the amplitudes Aa°\ Aa""^ 
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and Aa given by (|Tl|) with a = 1 for the plane layers 
(right frames) and a = 2 for the chain layers (left frames). 
For the orthorhombic system studied here all three of 
these harmonics belong to the same irreducible repre- 
sentation of the crystal point group except for 1712 = 
when the tetragonal Cu02 layer is decoupled from the or- 
thorhombic CuO layer. Recent cits measurements ||2^] 
show that the gap has a magnitude of ^ 20meV in the 
chains and ~ 30meV in the planes which would indicate 
that gi2 is large. 

For the first choice of intralayer interaction (upper 
frames), gn = 26.2, and there is no order parameter 
in the chains when there is no interlayer interaction (ie, 
.912 ~ 0) and the order parameter in the planes is pure d- 
wave. As the interlayer interaction is increased from zero, 
s-wave components appear in the planes and all three 
components appear in the chains. This "s + d mixing" is 
caused by the breaking of the tetragonal symmetry upon 
the introduction of the chains; there is no relative phase 
between the s- and d-wave components within either the 
planes or chains but there can be a relative phase be- 
tween the order parameter in the planes and chains. In 
the range of 512 explored here the c?-wave component in 
the plane remains dominant but for sufficiently strong in- 
terlayer interaction the isotropic s-wave component even- 
tually dominates |2^ (ie, the gap nodes disappear). For 
interaction parameters 312, 522} — {26.2,10,0} the 
critical temperature is lOOK and the maximum value of 
the gap in the Brillouin zone is 27.5meV in the planes 
and S.OmeV in the chains, while the maximum values on 
the Fermi surfaces are approximately 22meV and 7meV 
respectively. The ratio 2Z\niax/7c is 6.4 in the planes and 
1.9 in the chains. 

For the second choice of intralayer interaction (lower 
frames), gn — 9.18, there is no order parameter in either 
the chains or the planes when there is no interlayer inter- 
action (ie, gi2 = 0). As the interlayer interaction is in- 
creased d-wave and then s-wave components of the order 
parameter appear in both the planes and chains. Again, 
there is no relative phase between the s- and d-wave com- 
ponents within either the planes or chains but there can 
be a relative phase between the order parameter in the 
planes and chains. In the range of 1712 explored here, 
the d-wave component is dominant. At approximately 
912 = 15 the gap nodes no longer cross the Fermi sur- 
face in the chains; the feature at gi2 ^ 25 coincides with 
the gap nodes leaving the Brillouin zone and isotropic 
s-wave becoming dominant. For interaction parameters 
{5117512,522} = {9.18,20,0} the critical temperature is 
again lOOK and the maximum value of the gap in the Bril- 
louin zone is now 32.8meV in the planes and 20.1meV in 
the chains, while the maximum values on the Fermi sur- 
faces are approximately 27meV and 17meV respectively. 
The ratio 2Z\max/?c is 7.6 in the planes and 4.7 in the 
chains. 

Note that for 1712 > all of the s-wave components of 
the order parameters in both the planes and chains have 



the same relative sign and the d-wave components have 
opposite signs, while for 1712 < all of the relative signs 
are reversed but that the magnitudes of the components 
are insensitive to the sign of gi2 as noted after Eq. |^. 

In Fig. H we plot the density of states (dos) for the 
planes (left frames) and chains (right frames) calculated 
using the lowest three harmonics ( pi] ) of the solution to 
the BCS equations (^ with two sets of interaction pa- 
rameters. The dotted curves are the normal state DOS 
(nsdos) and the solid curves are the superconducting 
DOS (SCDOS). The insets show the Fermi surface (dashed 
curves) and gap nodes (solid curves) in the first Brillouin 
zone (with (7r,7r) at the center). The peak in the NS- 
DOS (dotted curves) is the van Hove singularity located 
at 2ta{2 — — 4:Ba ± Eq) which is at -62meV in the 
plane layers (left frames) and at 10 and -170meV in the 
chain layers (right frames). They are caused by the sad- 
dle points in the electron dispersions, ek.a, at (0, ±7r) and 
(±7r, 0). These van Hove singularities are shifted by the 
presence of the superconducting order parameter (solid 
curves) by an amount that depends upon the value of 
the order parameter at the saddle points; p^ ] in frame 
(c) these values are very different and the van Hove sin- 
gularity is split, in frame (a) these values are almost the 
same and no splitting is evident. 

An interesting feature is that the low energy behaviour 
{uj ^ 0) of the SCDOS is governed by the smallest local 
maxima of the gap on the Fermi surface when the gap 
nodes cross the Fermi surface and by the minima of the 
gap on the Fermi surface when there is no gap nodes 
which cross the Fermi surface. In Fig. ^ the magnitude 
of the gap along the Fermi surface is plotted. In Fig. IJa-c 
one can see that there are two different local maxima of 
the gap on the Fermi surface and these maxima are (to 
first order |^^) manifested as twin peaks in the SCDOS 
(Fig. ||a-c); these peaks are distinct from the van Hove 
singularities which are also present in the normal DOS 
(dotted curves in Fig. ^ and which are slightly shifted in 
the superconducting state. In frame (a) of Fig. ^ the 
local maxima of the gap on the Fermi surface are 16 and 
18meV; in (b) they are 1 and 7meV, and in (c) they are 
25 and 3meV. In (d) one can see that there are no gap 
nodes which cross the Fermi surface; the maximum and 
minimum value of the gap on the Fermi surface are 17 
and 4meV respectively. In Fig. ||d the finite gap in the 
SCDOS corresponds to the minimum of the gap on the 
Fermi surface and the peak to the maximum. 

The Josephson junction resistance-tunneling current 
product, RJ{T = 0), for a c-axis YBCO-Pb junction, 
given by ( [l^ ) with and appropriate for Pb, [H] are 
±0.25meV and ±2.2meV for the planes and 2.2meV and 
3.4meV for the chains for the two choices of gap made, 
in agreement with earlier calculations, jl^ The relative 
sign is due to the relative sign of the s-wave components 
(ie, the only part which contributes) of the order param- 
eters. The actual c-axis Josephson junction resistance- 
tunneling current product, RJ(T), for a junction made 
with untwinned YBCO would be some weighted average 
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of the plane and chain c-axis tunnehng currents depend- 
ing upon the relative abundance of chains and planes in 
the top layer of the YBCO. For a twinned sample with 
both twins equally abundant there would be zero net tun- 
neling current, although there is evidence that for single 
crystals of YBCO there can be up to a 5:1 ratio in the 
relative abundance of the two twin orientations. |^ We 
note that due to the different magnitudes of the order 
parameters in the two layers the model presented here 
is consistent with the observed tt shifts in corner junc- 
tions pO|-|3^ for both attractive and repulsive interlayer 
interactions, 512. 

In Fig. ^ we have plotted the magnetic penetration 
depth (left frames) and the Knight shift (right frames) 
calculated with the lowest three harmonics (^l]) of the 
solutions of the BCS equations (^) for the two choices of 
interaction parameters. In the penetration depth frames 
(left) the solid curve is for the x-direction (along the 
chains) and the dashed curve is for the y-direction (per- 
pendicular to the chains). The dotted curve is 1 — (T/Tc)^ 
and is plotted for comparison. The ratio \yy/Xxx a-t 
zero temperature is 1.37 for both interaction parameter 
choices since the zero temperature penetration depth is 
a normal state property (ie, the second term in Eq. 
does not contribute at zero temperature) . The zero tem- 
perature penetration depth is largely governed by the 
bandwidth (ie, 4ia(2 — e^)) - the larger the bandwidth 
the larger the zero temperature penetration depth. 

As pointed out above, the curvature of the penetra- 
tion depth curve, A~^^(r), is largely governed by the ratio 
2Z\max/Tc and is a straight line for the d-wave BCS value 
of 4.4. The presence of the chain layer and the interlayer 
interaction increases this ratio in the plane layer but it 
remains low in the chain layer due to the absence of an 
interaction in this layer. It is this lower value that makes 
Xyy{T) (along the chains) have upward curvature (solid 
curves). 

One can see (Fig. ^,c) that the in-plane penetra- 
tion depth perpendicular to the chains (dashed curve) 
closely resembles that observed experimentally in high 
quality crystals, 0j2^ and is largely determined by the 
presense of gap nodes crossing the Fermi surface, which 
cause the low temperature linear behaviour, and the ra- 
tio 2Z\niax/7c which, for values above ~ 4.4, make the 
curve of A~^^(r) have downward curvature. The pene- 
tration depth along the chains (solid curves), however, 
has an overall upward curvature due to the low values 
of 2Z\inax/7c in the chains. It is the component of the 
penetration depth due to the chains that makes the over- 
all \yy[T) have downward curvature. The component of 
the penetration depth due to the chains perpendicular to 
the chains does not contribute significantly to the overall 
penetration depth in this direction, )^yy , is due almost 
entirely to the Cu02 layers. For a single plane model 
|§ we would have {511,512,522} - {30,0,0}, 2A^,^/T, 
would be 4.5, and the penetration depth would closely 
resemble the straight line 1 - (T/Tc). Q 

The Knight shift (Fig. ||b,d), which is calculated in- 



dependently for the planes (solid curves) and chains 
(dashed curves), has a low temperature power law be- 
haviour when the gap nodes cross the Fermi surface 
(planes, both figures and chains in upper figure) and an 
exponential behaviour when the gap is finite over all the 
Fermi surface (chains, lower figure). When these quan- 
tities are measured the distinction between a power 
law and exponential behaviour rests upon the choice zero 
and so is not a reliable indicator of the presence of gap 
nodes on the Fermi surface. 



IV. CONCLUSION AND DISCUSSION 

We have derived a general expression for the Hamilto- 
nian in a multi-layer system and then made a simplifica- 
tion and have explicitly diagonalized the Hamiltonian. A 
set of two coupled bcs equations are then derived for this 
simplified system which is subsequently solved numeri- 
callyby a fft technique. This technique, unlike others, 
|p^ , P^ makes no assumptions about the functional form 
(and hence the symmetry) of the order parameter in ei- 
ther layer nor any relationship between the order param- 
eters in the different layers. 

Using the three lowest harmonics (^l]) of the solu- 
tions found for the coupled BCS equations (^) the mag- 
netic penetration depth (|l^), normal and superconduct- 
ing density of states ( |14|J15D , Knight shift ( |l3| ) and 
c-axis Josephson resistance-tunneling current products 
(|l6| ) were calculated. 

Our choice of electron dispersion relations were made 
so as to approximate the YBCO system in which there 
are layers consisting of Cu02 planes as well as layers 
which contain CuO chains. Our choice of interactions 
was made so that there is a pairing interaction in the 
Cu02 layer as well as an interlayer interaction, but no 
pairing interaction in the CuO layer. 

The solution of the BCS equations is predominately of 
a d-wave character, but because the tetragonal symme- 
try is broken by the presence of the chains there is some 
mixture of s-wave order parameter with no relative phase 
between the components, although the relative sign of 
the order parameter in the planes and chains may be 
±1. Further, we find that due to a symmetry in the set 
of coupled BCS equations derived, the sign of the inter- 
layer interaction affects only the relative sign of the order 
parameter in the two layers and not their absolute mag- 
nitudes (although some properties could be affected by 
this relative sign). We also find that any interlayer in- 
teraction strongly enhances the zero temperature order 
parameter, and hence the critical temperature. This is 
consistent with the observation that Tc is higher in ma- 
terials with multiple adjacent Cu02 layers. 

Furthermore, we find that the presence of gap nodes 
in only one of the layers is enough to produce a low tem- 
perature linear behaviour for the penetration depth, al- 
though if the minimum gap in the chains is too large the 
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A~^(T) and ^yy{T) curves can cross. Our calculation of 
the magnetic penetration depth gives a form that is simi- 
lar to that measured experimentally perpendicular to the 
chains, but not along the chains due to the small value of 
2Z\max/7c in the chains. This leads us to speculate that 
there may be intrinsic pairing in the CuO chains of the 
same order as in the Cu02 planes since the 2Z\inax/Tc ra- 
tios must both be large (6 to 7) for the low temperature 
slope of the A^^(T) curves to be approximately equal 
as is observed in experiments. This would tend to 

support the simple single band orthorhombic model pre- 
viously proposed by us. lTl|j20|] 

Our calculation of the superconducting density of 
states indicate that a surface probe may measure very 
different results depending upon whether the top layer 
is CuO chains or Cu02 planes. Depending upon the in- 
terlayer pairing strength the CuO chain layer may have 
a very narrow "c?-wave" type gap or a finite "isotropic 
s-wave" type gap. 

Our calculation of the c-axis Josephson resistance- 
tunneling current products, RJ{T = 0), for a YBCO-Pb 
junction for several choices of gap range from O.lSmeV 
to O.SOmeV for the planes and 2.35meV to 3.06meV for 
the chains, with possibly a relative sign between the 
chain and plane layers due to the relative sign of the 
s-wave components of the order parameters. The actual 
c-axis Josephson resistance-tunneling current product for 
a junction made with untwinncd YBCO would be some 
weighted average of the plane and chain results depend- 
ing upon the relative abundance of chains and planes in 
the top layer of the YBCO. For a twinned sample with 
both twins equally abundant there would be zero net tun- 
neling current although there is evidence that for single 
crystals of YBCO there can have up to a 5:1 ratio in the 
relative abundance of the two twin orientations. |29|1 



ACKNOWLEDGEMENTS 

Research supported in part by the Natural Sciences 
and Engineering Research Council of Canada (nserc) 
and by the Canadian Institute for Advanced Research 
(ciar). We would like to thank W. A. Atkinson and A. 
M. Westgate for discussions and insights. 



Electonic address: odonovan@mcmaster.ca 
[1] D.N. Basov, R. Liang, D.A. Bonn, W.N. Hardy, B. 

Dabrowski, M. Quijada, D.B. Tanner, J. P. Rice, D.M. 

Ginsberg and T. Timusk Phys. Rev. Lett. 74, 598 (1995). 
[2] A.J. Millis, H. Monien and D. Pines, Phys. Rev. B 42, 

167 (1990). 

[3] H. Chi and J.P Carbotte, Physics C 169, 55 (1990). 



[4] C. O'Donovan and J.P. Carbotte, Physica C 252 87 

(1995) . 

[5] D.A. Bonn, P. Dosanjh, R. Liang and W.N. Hardy, Phys. 

Rev. Lett. 68, 2390 (1992). 
[6] W.N. Hardy et al, Phys. Rev. Lett. 70, 3999 (1993). 
[7] D.A. Bonn, S. Kamal, K. Zhang, R. Liang, D.J. Baar, E. 

Kleit and W.N. Hardy, Phys. Rev. B 50, 4051 (1994). 
[8] N. Bulut and D.J. Scalapino, Phys. Rev. B 45, 2371 

(1992). 

[9] H. Ding, J.C. Campuzano et al., Phys. Rev. Lett. 50, 
1333 (1994). 

[10] H.L. Edwards, J.T. Markert, and A.L. de Lozanne, Phys. 

Rev. Lett. 69, 2967 (1992). 
[11] C. O'Donovan, D. Branch, J.P. Carbotte and J. Preston, 

Phys. Rev. B 51, 6588 (1995). 
[12] A.G. Sun et al., Phys. Rev. Lett. 72, 2267 (1994). 
[13] W.A. Atkinson and J.P. Carbotte, Phys. Rev. B 52, 

10601 (1995). 

[14] R.A. Klemm and S.H. Liu, Physics C 191, 383 (1992). 
[15] S.H. Liu and R.A. Klemm, Phys. Rev. B 45, 415 (1992). 
[16] R. Combescot and X. Leyranas, Phys. Rev. Lett. 75, 
3732 (1995). 

[17] S. Kettemann and K.B. Efetov, Phys. Rev. B 46, 8515 

(1992) . 

[18] T. Xiang and J.M. Wheatley, Phys. Rev. Lett. 76 134 

(1996) . 

[19] D.Z. Liu, K. Levin and J. Maly, Phys. Rev. B 51, 8680 
(1995). 

[20] C. O'Donovan and J.P. Carbotte, Phys. Rev. B 52, 4548 
(1995). 

[21] I. Maggio-Aprile, Ch. Renner, A. Erb, E. Walker and O. 

Fisher, Phys. Rev. Lett. 75 2754 (1995). 
[22] T. Jacobs, S. Sridhar, Qiang Li, CD. Gu and N. 

Koshizuka, Phys. Rev. Lett. 75, 4516 (1995). 
[23] Mechanisms of Conventional and High Tc Superconduc- 
tivity, V.Z. Kresin, H. Morawitz and S.A. Wolf, Oxford 

University Press (1993). 
[24] V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 

486 (1963); 11, 104 (1964). 
[25] K. Gofron, J.C. Campuzano, A. A. Abrikosov, M. Lin- 

droos, A. Bansil, H. Ding, D. Koelling, and B. Dabrowski, 

Phys. Rev. Lett. 73, 3302 (1994). 
[26] H.L. Edwards, D.J. Derro, A.L. Barr, J.T. Markert, and 

A.L. de Lozanne, Phys. Rev. Lett. 75, 1387 (1992). 
[27] A.L Liechtenstein, I.I. Mazin and O.K. Anderson, Phys. 

Rev. Lett. 74, 2303 (1995). 
[28] C. Zhou and H.J. Schuiz, Phys. Rev. B 45, 7397 (1992). 
[29] B. Gaulin, private communication. 

[30] D.A. Wollman, D.J. Van Harlingen, W.C. Lee, D.M. 
Ginsberg and A.J. Leggett, Phys. Rev. Lett. 71, 2134 

(1993) . 

[31] D.A. Brawner and H.R. Ott, Phys. Rev. B 50, 6530 

(1994) . 

[32] C.C. Tsuei et al, Phys. Rev. Lett. 73, 593 (1994). 
[33] A. Mathai, Y. Gim, R.C. Black, A. Amar and E.G. WeU- 
stood, Phys. Rev. Lett. 74, 4523 (1995). 



7 




1 
1 

/ . 

/ / 






1 
1 
\ 

^\ 1 - 
\-\ - 


// 

1/ 
w 

- I-! 






W - 

w 

I 


11 

. W 

K 






l\ 
II 

// " 


\ \ 
1 ^- 

1 

1 






/-/ - 

y 1 
1 

1 , " 



FIG. 1. Model of YBCO Fermi surfaces for chains (long 
dashed curve) and planes (closed short dashed curve) in the 
first Brillouin zone. The {n,n) point is at the center of 
the figure. For the chains the parameters {to , , , Ma } 
in Eq. |l^ are {—50, —0.9, 0, 1.2} and for the planes they are 
{100,0,0.45,0.51}. 
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FIG. 2. Calculation of the zero temperature order param- 
eters as a function of the interlayer interaction, (;i2, for two 
fixed values of the interlayer interaction, gn, (upper and lower 
frames) presented for the planes (right frames) and chains 
(left frames) separately. In all frames the solid curve is the 
d-wave component of the order parameter, the short dashed 
curve is the extended s-wave component and the long dashed 
curve is the isotropic s-wave component. In the upper frames 
gii = 26.2 and for g^i = 10, Tc = IQQK. At gu = the 
order parameter is zero in the chains and is pure d-wave in 
the planes. As the interlayer interaction is increased the or- 
der parameter becomes present in the chains and there is a 
mixing of s-wave components. In the lower frames = 9.18 
and for g\2 = 20, Tc = IQQK. At g\2 ~ the order parameter 
is zero in both the chains and the planes. As the interlayer 
interaction is increased the order parameter becomes present 
in both the chains and planes and there is a mixing of s-wave 
components with the isotropic s-wave component eventually 
becoming dominant. The feature at gi2 ~ 25 occurs when 
the gap node leaves the Brillouin zone. As discussed in the 
text there is a gi2 <-> —gi2 symmetry. Both gn and gi2 are 
in units of ti; (722, the coupling in the chains, is set equal to 
zero. 
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FIG. 3. Calculation of the superconducting (solid curves) 
and normal (dotted curves) density of states (dos) for two 
sets of interaction parameters (upper and lower frames) 
presented separately for the planes (left frames) and chains 
(right frames). Some experiments are surface probes and may 
probe either the planes or chains independently. The in- 
sets show the Fermi surface (dashed curve) and gap nodes 
(solid curve) in the planes and chains for the two dif- 
ferent parameter choices in the first Brillouin zone with 
(tt, 7r) at the center. In the upper frames, (a) and (b), 
we have chosen {(/n, (;i2, (?22} = {26.2,10,0} which gives 
Tc = lOOif. In the lower frames, (c) and (d), we have cho- 
sen {511,312,922} = {9.18, 20,0} which also gives Tc = lOOif. 
Note that for the second parameter choice the gap nodes do 
not cross the Fermi surface in the chains (frame (d), inset) 
and that the DOS is gapped. The c-axis Josephson resis- 
tance-tunneling current product for a YBCO-Pb junction for 
a pure d-wave order parameter is zero due to the equal parts 
of the order parameter with opposite signs. Here this is not 
the case (insets) and the c-axis Josephson resistance-tunneling 
currents product for a YBCO-Pb junction are (a) ±0.18meV, 
(b) 2.35meV, (c) ±0.25meV, and (d) 2.17meV. The relative 
sign is due to the relative sign of the s-wave components of 
the order parameter (ie, the only part which contributes). 
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FIG. 4. The magnitude of the gap on the Fermi surface 
as a function of angle for the four cases of Figs. ^ and ^ The 
angle 9 is measured from the center or (tt, tt) point of the 
Brillouin zone with the y-axis (ie, the vertical in the insets of 
Fig. 1^ corresponding to S = 0. Frames (b) and (d) do not 
span all angles due to the Fermi surface not being closed in 
the chain layer. For the first choice of interaction parameters, 
{311. 912,522} = {26.2,10,0}, the ratio 2A^^^(^s) /T^, where 
^max(FS) is the maximum value of the gap on the Fermi sur- 
face, is 4.3 and 1.6 for the planes and chain respectively; for 
the second, {(711,912,922} = {9.18,20,0}, they are 5.7 and 
3.8. 
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FIG. 5. Calculations of the magnetic penetration depth 
(left frames) and the Knight shift (right frames) for the two 
sets of interaction parameters. Frames (a) and (c) show the 
magnetic penetration depth along (solid curve) and perpen- 
dicular (dashed curve) to the CuO chains. The dotted curve is 
1 - {T/T^f and is shown for comparison. The chains, due to 
their Fermi surface, do not contribute appreciably to the pen- 
etration depth perpendicular to the chains (dashed curves). 
The ratio {\yy / \xx)^ is 1.37 for both sets of parameters since 
this is a normal state property. Frames (b) and (d) show the 
Knight shift in the planes (solid curves) and chains (dashed 
curves); due to the crossing of the gap nodes and the Fermi 
surface in the chains in (b) the Knight shift in the chain is a 
power law at low temperature, while due to the finite gap in 
(d) the Knight shift in the chain is exponential at low tem- 
perature. 
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